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Abstract 

An infinite family of classical superintegrable Hamiltonians denned on the N- 
dimensional spherical, Euclidean and hyperbolic spaces are shown to have a com- 
mon set of (2N — 3) functionally independent constants of the motion. Among 
them, two different subsets of N integrals in involution (including the Hamil- 
tonian) can always be explicitly identified. As particular cases, we recover in a 
straightforward way most of the superintegrability properties of the Smorodinsky- 
Winternitz and generalized Kepler-Coulomb systems on spaces of constant curva- 
ture and we introduce as well new classes of (quasi-maximally) superintegrable 
potentials on these spaces. Results here presented are a consequence of the 
sZ(2,R) Poisson coalgebra symmetry of all the Hamiltonians, together with an 
appropriate use of the phase spaces associated to Poincare and Beltrami coordi- 
nates. 
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1 Introduction 



An JV-dimensional (iVD) completely integrable Hamiltonian is called maximally 
superintegrable (MS) if there exists a set of (2JV — 2) globally defined functionally 
independent constants of the motion that Poisson-commute with H^ N \ Among these 
constants, at least, two different subsets of (JV— 1) constants in involution can be found; 
these are connected with the fact that the system must possess separable solutions to 
the Hamilton- Jacobi equation in, at least, two different coordinate systems. From the 
dynamical point of view, all bounded motions of MS systems are closed and strictly 
periodic. In a similar way, if has (2iV — 3) integrals with the abovementioned 
properties the system will be called quasi-maximally superintegrable (QMS). Obviously, 
all MS systems are QMS ones. 

The search of explicit MS (or QMS) systems is usually afforded by fixing the number 
N of degrees of freedom together with some assumptions concerning both the functional 
dependence of H^ N ' and of the integrals of the motion. In the case of natural systems 
with H^ N > = T + U, the configuration space of the system (a space with a given 
curvature, for instance) is chosen through a fixed expression for the kinetic energy T, 
and MS potentials U are investigated by imposing the existence of a suitable number 
of (unknown) integrals of the motion, that are usually assumed to be quadratic in the 
momenta and have also to be found simultaneously. 

However, it turns out that although for a given space and dimension there exists 
a certain number of superintegrable potentials, only very few of them admit arbitrary 
iVD generalizations. For example, in the JVD Euclidean space the only two known 
MS Hamiltonians with integrals quadratic in the momenta are the oscillator potential 
with JV centrifugal terms (the so called Smorodinsky-Winternitz system PUI2]) an d the 
generalized Kepler-Coulomb potential with (JV — 1) centrifugal terms jSJ. 

The aim of this work is to face this problem from a quite different viewpoint, which 
is based on introducing a common si (2, M) <8> • • <8> sl(2, R) symmetry framework for a 
large family of JVD QMS systems with integrals quadratic in the momenta. Surprisingly 
enough, all the systems introduced in this way will share, by construction, the same set 
of (2JV — 3) functionally independent integrals of the motion, which will be explicitly 
given. Among them, only one constant of the motion C^ N ' will depend on the JV 
degrees of freedom and the remaining ones will be functions of a decreasing number of 
positions and momenta. 

Therefore, in this framework the Hamiltonian will be the only function that 
characterizes each individual QMS system. In the case of natural systems, the kinetic 
term T will give us the information concerning the configuration space, and the po- 
tential term U will fully identify the interactions. We will show that, by following this 
approach, QMS systems on JVD spaces of constant curvature can be easily understood 
and constructively studied. In this respect, the geometric interpretation of the canon- 
ical coordinates and momenta on each space turns out to be essential, and will be 
induced by the Hamiltonian transcription of the Poincare and Beltrami coordinates. 

Moreover, we stress that for some specific choices of the systems here presented 
are indeed MS ones, and we will recover as particular cases the JVD Smorodinsky- 
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Winternitz and generalized Kepler-Coulomb systems, that can now be understood as 
distinguished cases of an infinite family of systems with the same underlying symmetry. 
For both cases the additional independent integral of the motion, which is not provided 
by the underlying symmetry, can be found by other methods and will be explicitly 
given. 

The paper is organized as follows. In the next Section we present the main result 
that defines the infinite family of QMS systems together with the explicit form of their 
'universal' (2N — 3) integrals of the motion. A sketch of the proof of this theorem is 
given, which is based on the symmetry of our systems under an si (2, K.)®. . . N ^(g)sl(2, R) 
algebra (or, alternatively, in terms of an s/(2,M) Poisson coalgebra [HE]). Afterwards, 
many different examples of Hamiltonian systems belonging to this family (and thus 
sharing the same set of universal integrals) are presented. In Section 3 we construct 
QMS systems on the ND Euclidean space by using a kinetic energy T given in 
terms of the usual momenta conjugated to Cartesian coordinates, that can be naturally 
identified with the canonical coordinates q. Some of these systems were already given 
in El E| but others are here identified within the coalgebra symmetry framework 
for the first time. Section 4 is devoted to QMS systems living on the ND spherical 

and hyperbolic M N spaces, for which the coordinates q are identified either with 
Poincare or with Beltrami coordinates coming, respectively, from stereographic and 
central projections on M. + . Finally, some remarks and open problems are mentioned. 

2 An infinite family of QMS systems 

The main result of this paper can be stated as follows. 

Theorem 1. Let {q, p} = {(qi, . . . , gjv), (pi, ■ ■ ■ ,Pn)} be N pairs of canonical vari- 
ables. The ND Hamiltonian 

#W=ft(q 2 ,p 2 ,q-p) (1) 

with 7i any smooth function and 

N N / h \ N h N 

q 2 = E«? p 2 = £ P? + ? q-P=E*fl. (2) 

i=l 1=1 v q% ' 1=1 % i=l 

where bi are arbitrary real parameters, is quasi-maximally superintegrable. The (2N—3) 
functionally independent and 'universal' integrals of the motion are 

c {m) = £ [te - qjPt ) 2 + + 6i|j }+:>> 

N ( / 2 2\ ^ N 

C(m) = \ " + ( bA + & A ) | + b * ( 3 ) 

N-m+l<i<j ^ V % 1j J ) i=N _ m+1 

where m — 2, . . . , N. Moreover, the sets of N functions {H^ N \ C^™'} and {H^ N \ C( m )} 
(m = 2, . . . , N) are in involution. 
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Some remarks are in order: 



• Note that C^ N > = C(jv), so that the number of different integrals is (2N — 3). 

• The underlying symmetry can be interpreted as a generalization of the 'spher- 
ical' one, since H^ N ' = 7i(q 2 ,p 2 ,q ■ p) is recovered when all hi = 0. In this 
case, the constants of the motion @ are just sums of squares of certain angular 
momentum components kj = qiPj — qjPi. In particular, C^ m ' = YlT<i<j^ij anc ^ 
C — I 2 

Km) A^N-m+l<i<j 

• As we shall see in the sequel, the additional centrifugal terms h/q 2 in come 
from the freedom to add such a term in the corresponding symplectic realization of 
an si(2,R) Poisson algebra. These terms will be essential in order to understand 
several known results concerning QMS systems on spaces of constant curvature. 

Proof. This relies on the fact that, for any choice of the function 7i, the Hamiltonian 
H( N ) has an si(2,R) Poisson coalgebra symmetry [S] (under a certain symplectic re- 
alization) and, as a consequence, the integrals of the motion come from the left and 
right m-coproducts of the Casimir function for sZ(2, R) [HIE!- Let us summarize the 
main steps of this construction. 

We recall that the sZ(2,R) Poisson coalgebra is given by the following Lie-Poisson 
brackets, coproduct, Casimir and (primitive) coproduct A: 

{J 3 ,J + } = 2J + {J 3 ,J_} = -2J_ {J_,J + } = 4J 3 (4) 

A(J { ) = Ji ® 1 + 1 ® Ji / = +,-, 3 (5) 
C = J_J + -Jl (6) 
A one-particle symplectic realization of (J3J) is given by 

J- = q\ J+=p\ + b 1 /ql Js = qiPi (7) 

where b\ is a real parameter that labels the representation through C = b\. The 
coalgebra approach [S] provides the corresponding A^-particle symplectic realization 
through the iV-sites coproduct of © living on sl(2, R) ® . . ® sl(2, R) 0: 

jv n , , v N b N 

J- = ^q 2 = q 2 J+ = ( V\ + ~\ ) = P 2 + ~k Js = S QiPl = q ' P ( 8 ) 
i=l i=i \ ^ ' i=l ^ i=i 

where 6, are iV arbitrary real parameters. This means that the A^-particle generators 
(JHJ) fulfil the commutation rules (J3J) with respect to the canonical Poisson bracket 

N 



f-f\dqidpi dqidpij 



As a consequence of the coalgebra approach, these generators Poisson commute 
with the (2N — 3) functions Q given by the sets C^™' and C( m ), which are obtained 
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from the 'left' and 'right' m-th coproducts of the Casimir © with m = 2, 3, . . . , N 
(see [7j for details). If we label the N sites on sl(2,R) <g> sl(2,R) ® . . . g> sZ(2,R) by 
1 <£) 2 (g> . . . ® N, the 'left' Casimir C^ m ^ is defined on the sites 1 (g> 2 (g> . . . (g> m, while 
the 'right' one C( m ) is defined onm®...8iV- 1 ® iV. Moreover, it is straightforward 
to prove that the 2N - 2 functions {C (2 \ . . . , = C {N)} C {N _ 1)} . . . C {2) ,H} are 
functionally independent (assuming that 7i is not a function of C) and the coalgebra 
symmetry ensures that each of the subsets {C^, . . . , C^ N \H} or {C( 2 ), . . . , Cm), T~L} are 
formed by N functions in involution [13 [7|. 

Therefore, any arbitrary function 7i defined on the iV-particle symplectic realization 
of sl(2, R) © is of the form that is, 

^)=^(J_,J + ,J 3 )=H^q 2 ,p 2 + g|,q-pj 

and defines a QMS Hamiltonian system. 

Notice that when N = 2, the generic function 7i determines an integrable Hamil- 
tonian as there is a single constant of the motion = C(2). On the contrary, when 
N = 3 the Hamiltonian 7^ (0) can be called minimally or 'weak' superintegrable (by 
following the terminology of [3] for this specific dimension) since it is endowed with 
three integrals {C^ 2 \ C(2)}- However, we remark that for arbitrary iV there is a 
single constant of the motion left to assure maximal superintegrability. In this respect, 
we stress that some specific choices of 7i comprise maximally superintegrable systems 
as well, but the remaining integral does not come from the coalgebra symmetry and 
has to be deduced by making use of alternative procedures. 



3 QMS systems on the Euclidean space 

It is immediate to realize that the iVD kinetic energy on directly arises through the 
generator J + (jHJ) as it includes the usual momenta p 2 . Then if we add some smooth 
function J~(J-) as the potential term we can construct different QMS systems. Some 
possibilities are the following. 

• The ND Evans system: 

1 1 N h 

1=1 1 

where hi = bi/m. This describes a particle of mass m under the action of a central 
potential ^"(q 2 ) with iV centrifugal terms associated with the 6j's. The iV = 3 case was 
introduced in j^j. 

• The Smorodinsky-Winternitz system [TJ |2l El El EE UD] • By choosing J 7 (J J) = uj 2 J_ in 
the above Hamiltonian, we recover a system formed by an isotropic harmonic oscillator 
of mass m with angular frecuency uj together with iV centrifugal barriers: 

1 1 - h- 

n = ^-J + + ^ = J- p 2 + ^ 2 q 2 + J2 A- (io) 

2m 2m ' 2qf 

i=l 1 
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This Hamiltonian is known to be maximally superintegrable. A set of iV additional 
constants of the motion (which do not come from the coalgebra symmetry) is 



Ti = pf + 2muj 2 qt + mbi/qf i = l,...,N. (11) 

Each constant X, is functionally independent with respect to both the set @ and TC, 
thus it can be taken as the 'lost' integral. 

• A Gamier-type system |11| IT2*]. If we set T{JJ) = c<j 2 J_ + <5J 2 , with real parameter 5, 
we find a degenerate Gamier Hamiltonian (with = u 2 ), which besides the previous 
harmonic term also comprises quartic oscillators: 

i i N h 

H - i J + + " 2J - +SJ --^ 2+ +s « ,+ ^w < 12 > 

i=l 1 

The generalization to even-order nonlinear oscillators is straightforward, since the 
generic system 

n = i J+ + " 2J - + £* J - +1 = ip 2 +-V +£ + < 13 > 

k=l k=l i=l y * 

is a QMS one for any choice of the 5k parameters. 

— 1/2 

• A generalized Kepler-Coulomb system jo]. The choice F{J~) = —kJ_ , with real 
constant k, gives rise to the superposition of the Kepler-Coulomb potential with N 
centrifugal barriers: 

n = —j + - u- i/2 = — P 2 -^+y\ (14) 



2m 2m ^£ f-f 2q, 

Such a Hamiltonian is known to be maximally superintegrable under the condition 
that, at least, one of the centrifugal terms vanishes. In particular, if a single 6j = 0, 
an additional constant of the motion (and functionally independent with respect to 
and JHJ) arises. This additional integral is found to be 

N N 

1i i Qi 



Ct = YVftft - ftp,) +to^=-mr^. (15) 



Likewise, if another bj = (j ^ i), the function £j can be proven to be also a constant 
of the motion, and so on. When all the fej's vanish the system reduces to the proper 
Kepler-Coulomb potential and the N additional constants of the motion Li are the 
components of the Laplace-Runge-Lenz iV-vector on E . 

• Stationary electromagnetic fields. Certain momenta-dependent potentials can also be 
obtained within this framework through an additional term depending on the generator 
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Js- Let us consider 



n = ^-j + --j 3 g( j_) + ej^( j_) 

2m m 

= f- P 2 - - (P " q)^(q 2 ) + e-^(q 2 ) + £ £ (16) 

i=i ^ 

where e is a real parameter and Q{J~) is a smooth function. When iV = 3, the above 
Hamiltonian describes a particle with mass m and charge e that moves on E 3 under 
the action of an electromagnetic field, that is, TC = j^(p — eA) 2 + eip, where the 
(time-independent) vector A and scalar ip potentials are given by 

V>(q) = ^(q 2 ) - j- q 2 £ 2 (q 2 ) + £ ^ A(q) = q£(q 2 ). 

" i=l * 

The electric E = — V"0 and magnetic H = V x A fields turn out to be 

B=f± < r + ^ <rt( r-2r)q+±&*»^ H = 



m m / e^Q'i Q'i % 

where Q' and J 7 ' are the derivatives with respect to their variable q 2 . Recall that 2D 
integrable electromagnetic Hamiltonians have been studied in [T3*l ITU ITS] . We stress 
that this kind of construction can also be applied to obtain ND QMS Fokker-Planck 
Hamiltonians (see and references therein). 

• Systems with coordinate- dependent mass. If we multiply the former kinetic term J + 
in the above Hamiltonians by an arbitrary smooth function Ai(J_), we obtain systems 
with a variable mass depending on q 2 such as, for instance, 

1 1 N h 

n -rnm^ +nj -^w/ +n ^%^wr (17) 

We stress that the expression (|T7jl can be alternatively interpreted as a QMS system 
on 'some' space which would be determined through the metric coming from the kinetic 
energy J + / 'Ai(J-). In particular, in the next Section we will show how certain specific 
Hamiltonians of the type (|T7|) will define superintegrable systems on the sphere E> N and 
the hyperbolic space M N . 



4 QMS systems on the sphere and the hyperbolic 
space 

Recall that for a constant sectional curvature k, both the ND sphere (k > 0) and 
the iVD hyperbolic space H w (k < 0) can be embedded in an ambient linear space M Ar+1 
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with Weiertrass coordinates (xo,x) = (xq,xi, . . . ,xn) fulfilling the 'sphere' constraint 
S: x\ + kx 2 = 1. The metric on the proper iVD spaces is given by [THJ [T7j: 

ds 2 = - (dxl + Kdx 2 ) (18) 

K £ 

where dx 2 = ^ i=1 da; 2 . The N + 1 ambient coordinates can be parametrized in terms 
of N intrinsic coordinates in different ways. 

Firstly, we consider the stereographic projection JH] from the ambient coordinates 
(x ,x) 6 E to the Poincare ones yet" with pole (-1, 0) G R N+1 . Hence (-1, 0) + 
A (1, y) G S and we obtain that 

A = — 2 x = A - 1 = — ^ x = Ay = — ^ 19 

1 + Ky z 1 + ny z 1 + ny z 

Secondly, we apply the central projection from the Weiertrass coordinates to the Bel- 
trami ones z G M. N with pole (0, 0) G R N+1 . Thus (0, 0) + // (1, z) G £ and we find 

/i = - 1 = x = /j x = /iz= - Z = . (20) 



KZ 



Then the metric ()18j) in Poincare and Beltrami coordinates turns out to be 

j 2 _ a dy 2 (1 + /€Z 2 )dz 2 - k(z ■ dz) 2 

(l + /ty 2 ) 2 (1 + kz 2 ) 2 ' 1 J 



Next if we write the metrics as Poincare T p and Beltrami T B free Lagrangians 

P = my 2 t b = m / (1 + /tz 2 )z 2 - /c(z • z) 2 

2(l + /ty 2 ) 2 2 V (1 + kz 2 ) 2 

the momenta p y and p z conjugated to y and z can be deduced 



my / (1 + kz 2 )z — k(z ■ z)z\ 

~ m ~ T, ..A ■ (23) 



^ y (1 + ^y 2 ) 2 cz V (l + /tz 2 ) 2 

Therefore, the previous expressions show that the ND kinetic energy in both 
Poincare T p and Beltrami T B phase spaces can be written as the following functions 
of the symplectic realization (JBJ) of the generators of si(2, R): 

rP = i (1 + Kj - )2j+ = i( 1 + Kq2 > v 

Am Am (24) 

r B = ^ (i + kJ-) (j+ + « Jf ) = 2^(1 + ^q 2 ) (p 2 + «(q ■ p) 2 ) 

provided that all fej = 0, (q, p) = (y, p y ) for T p and (q, p) = (z, p z ) for T B . When the 
fej's are taken as arbitrary parameters, the N additional curved 'centrifugal' potentials 
arising in (J24|) are written in ambient coordinates in a very simple form: 

N h N h N h N h 

^ 2my 2 V y ' x 2 2mz 2 V > ^ 1x 2 K ' 

i=l y * i=l * »=1 * i=l 1 



tan 2 (v^ r) = -5 = _ = z 2 . (26) 



It is worth remarking that (|23j) is actually a direct consequence of the symplectic real- 
ization of sl(2, M) (JSJ) with 6j 7^ 0. This, in turn, explains why the famous 'centrifugal' 
terms do appear everywhere in the classifications of superintegrable systems on the 
three classical Riemannian spaces UHl 1201 EEj • 

Recall as well that the radial (geodesic polar) distance r from an arbitrary point to 
the origin in and H w along the geodesic joining both points is expressed in terms 
of Weiertrass ^01 HZj , Poincare and Beltrami coordinates as 

1 x 2 4y 2 

^ tan2( ^ r) = ^ = (T^ 

Note that in the constant curvature analogues of the oscillator and Kepler-Coulomb 
problems the Euclidean radial distance is just replaced by the function tan(-y/rer). 

Therefore, by taking into account (|24 )) -([26 )) we finally present some QMS and MS 
Poincare 7i p and Beltrami 7i B Hamiltonians that are constructed by adding some 
suitable functions depending on J_ to (J24|) and by considering arbitrary centrifugal 
terms Vs. These systems are the curved counterpart of the Euclidean systems (jHJ)- 
(|T3j) which simultaneously cover (re > 0), H N (re < 0), and K N (k, = 0; in this flat 
limit r 2 = x 2 = 4y 2 = z 2 ). We stress that, again, all of them share the same set of 
constants of the motion (JHJ), although in all these curved cases the geometric meaning 
of the canonical coordinates and momenta is completely different to the Euclidean one, 
as explained above. To make this geometrical interpretation more explicit, for each 
system we shall give both the Poincare and Beltrami versions of each Hamiltonian. 



A curved Evans system. It would be given by 

^P_ T P | T ( 4J_ Wl + .q 2 )V / 4q 2 \ *2 bl 

^ Wl-reJ_W 2m W 1 - req 2 ) 2 / ^-f 



.v V =1 ' (27) 



<H* = T B + T ( J.) = ±{1 + .q 2 ) (p 2 + re(q ■ p) 2 ) + T (q 2 ) + £ ^. 

i=l % 

• TTie curved Smorodinsky-Winternitz system [TBJ[T7|. Such a system is just the Higgs 
oscillator [221 122] (that arises as either polynomial or rational potential in these coor- 
dinates) plus the corresponding centrifugal terms: 

n r _ t p , 4^ 2 J. _ (1 + Kg 2 ) 2 p 2 | 4^ 2 q 2 | " 2b t 
(1 - kJ_) 2 2m (1 - Kq 2 ) 2 ^ x\ 



i=l 

N 



H» = T* + ^J. = ±{1 + .q 2 ) (p 2 + re(q ■ p) 2 ) + .V + Y: 2j , • 

~" ' i=l 1 



{21 



This Hamiltonian is maximally superintegrable and the remaining constant of the mo- 
tion can be taken from any of the following N functions (to be compared with 

^■v i t , 2\ r. / \ \2 8mu 2 qf ~ (1 - /tq 2 ) 2 

If = ( Pi (l - K q 2 + 2re q ■ p ft ) + ^ + mfe, ^ ^ 

(1 - nq 2 ) 2 qf (29) 

2? = (Pi + K (q • P)<?i) 2 + 2ma; 2 g 2 + mbi/qj i = 1, . . . , N. 
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• A curved Gamier-type system. By following the same prescription, QMS quartic 
curved oscillators can be defined as 

H p = T - + ,. 4 " V - + ^ 



A-K.J-) 2 (1-/«J_) 4 

(l + /tq 2 ) 2 p 2 Au\ 2 16£q 4 A 2h 

2m + (1 - /tq 2 ) 2 + (1 - K q 2 ) 4 + ^ xf ( 30 ) 

l + Kq 2 )(p 2 + K(q-p) 2 ) , 2 , , , 4 ( A bi 



2m ^ 2x 2 

i=l 1 

The expressions for the curved QMS analogues of the higher-order nonlinear oscillators 
f]"3|) are straightforward. 

• A curved generalized Kepler-Coulomb system I2DI EH I2H E3 The curved 
Kepler-Coulomb potential with N centrifugal terms corresponds to 



4J_ \ - 1/2 (1 + Kq 2 f p 2 , (1 - /tq 2 ) A 2b t 



Wl- K J_ 2 i 2m 2 A AF 4f x 2 



q 

1 i=i 

i k b- 



(31) 



^ B = T B - kJZ 1/2 = + .q 2 ) (p 2 + «(q ■ p) 2 ) -^= + £^2- 



This is again a MS system provided that, at least, one bi = 0. In this case the remaining 
constant of the motion reads (compare with the Euclidean version ([15)1 ): 

A kmq { A ~ &(1 - Kq 2 ) 



^ 2 



n kmq- N 



«=1 V "1 i=l;/^j 11 

If another 6j = 0, then £j' B is also a new constant of the motion. In this way the 
proper curved Kepler-Coulomb system [27] (with all the 6;'s equal to zero) is obtained, 
and in that case (|32JI are the N components of the Laplace-Runge-Lenz vector on 
(k > 0) and M N (k < 0). 

More details on the latter Hamiltonians and their generalization to Lorentzian met- 
rics will be given elsewhere. On the other hand, we remark that the definition of QMS 
systems on spaces of variable curvature can be achieved either by considering more gen- 
eral kinetic energy terms or by making use of coalgebra deformations, that have been 
already shown to underly the superintegrability on some variable curvature (Rieman- 
nian and relativistic) spaces [2E1I2H]- Finally, the study of several interesting classes of 
non-natural Hamiltonians included in ((TJ) are also worthy to be addressed in the future. 
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